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In this paper we measure the angular power spectra Ce of three high-redshift large-scale structure 
probes: the radio sources from the NRAO VLA Sky Survey (NVSS), the quasar catalogue of Sloan 
Digital Sky Survey Release Six (SDSS DR6 QSOs) and the MegaZ-LRG (DR7), the final SDSS II 
buminous Red Galaxy (bRG) photometric redshift survey. We perform a global analysis of the 
constraints on the amplitude of primordial non-Gaussianity from these angular power spectra, as 
well as from their cross-correlation power spectra with the cosmic microwave background (CMB) 
temperature map. In particular, we include non-Gaussianity of the type arising from single-field 
slow roll, multifields, curvaton (local type), and those which effects on the halo clustering can be 
described by the equilateral template (related to higher-order derivative type non-Gaussianity) and 
by the enfolded template (related to modified initial state or higher-derivative interactions). When 
combining all data sets, we obtain limits of /nl = 48 ± 20, /nl = 50 ± 265 and /nl = 183 ± 95 
at 68% confidence level for local, equilateral and enfolded templates, respectively. Furthermore, we 
explore the constraint on the cubic correction (?nl0 3 on the bias of dark matter haloes and obtain 
a limit of -1.2 x 10 5 < #nl < 11.3 x 10 s at 95% confidence level. 



I. INTRODUCTION 

The physical mechanisms responsible for the generation of primordial perturbations seeding present-day large-scale 
structure, may leave their imprint in the form of small deviations from a Gaussian distribution of the primordial 
perturbations. Searches for primordial non-Gaussianity can thereby provide key information on the origin and evo- 
lution of cosmological structures (e.g., ref. pj and references therein). Although the standard single-field, slow- roll, 
canonical kinetic energy and adiabatic vacuum state inflation generates very small non-Gaussianity, any inflationary 
model that deviates from this may entail a larger level of it (refs. @, Q an d references therein). 

Deviations from Gaussian initial conditions arejaften taken to be of the so-called local type and are parameterized 
by the constant dimensionless parameter /nl 



$ = </> + /nl(0 2 -(0 2 )) , (1) 

where $ denotes Bardeen's gauge-invariant potential (evaluated deep in the matter era) and </> is a Gaussian random 
field 1 . 

Recently, a method for constraining non-Gaussianity from large-scale structure surveys has been proposed fol. [Toj. 
which exploits the fact that the clustering of dark matter halos -where galaxies form- on large scales is modified in 
a scale-dependent way by the presence of even small amount of non-Gaussianity. In particular, a non-Gaussianity 
described by cq. (JlJ, introduces a scale-dependent boost (for /nl > and a suppression for /nl < 0) of the halo 
power spectrum proportional to 1/fc 2 on large scales (k < 0.03 h/Mpc), which evolves roughly as (1 + z). Lar ge-S cale 
Structure (LSS) surveys covering large volumes are needed to access the scales where the signal arises (e.g., ref. [l2| and 
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references therein). The possibility of using high redshift data to constrain non-Gaussianity has also been addressed 
by means of hydrodynamic simulations of the Intergalactic Medium [l3j . Among the many currently-available tracers 
of the LSS, the radio sources from the NRAO VLA Sky Survey (NVSS) Q and the quasar catalogue of the Sloan 
Digital Sky Survey Release Six (SDSS DR6 QSOs) [15| are particularly interesting since they span large volumes 
extending out to substantial redshifts [l6l ]. Indeed these source samples were shown to provide tight constraints on 
primordial non-Gaussianity by the cross-correlation measurements between CMB temperature fluctuations and the 
LSS number density which could be used to detect the late Integrated Sachs- Wolfe (ISW) [l7| effect [I8I - I2TI ]. 

In recent papers 0, [2l| it was shown that the observed NVSS and SDSS DR6 QSOs Auto-Correlation Function 
(ACF) and Cross- Correlation Function (CCF) hint at a positive value of /nl at more than 2 er. This is because the 
ACF is found to be still positive on angular scales 9 > 4°, which, for the median source redshift (z m — 1), correspond 
to linear scales where the correlation function should be negative if the density fluctuation field is Gaussian. A positive 
/nl adds power on large angular scales, accounting for the observed ACF. 

In this paper, we will revisit the constraints on the primordial non-Gaussianity from these two LSS tracers using 
the power spectrum Ce, instead of the correlation function w(8). While it is true that the angular correlation function 
and the angular power spectrum arc spherical harmonic pairs and thus in principle they should enclose the same 
physical information, each of the two measurements is affected by systematic errors in different ways. For example 
the correlation function is less affected by the survey mask than the power spectrum but error estimation for the power 
spectrum is much easier then for the correlation function. A direct comparison of parameter constraints obtained via 
the correlation function and the power spectrum analyses offer therefore a way to quantify possible systematic effects. 
We will also use the recently updated SDSS II LRG photometric redshift survey, MegaZ-LRG (DR7) [22| to constrain 
the primordial non-Gaussianity. Ref. [l8| used the LRG angular power spectrum for a sample close to the SDSS 
DR3 to constrain /nl- Here we revisit the analysis using an improved bias model, improved LRG catalog, updated 
knowledge of the sources redshift distribution, up-to-date CMB maps and complementary analysis methods. 

Furthermore, we explore the constraints on physically motivated primordial non-Gaussianity shapes that are dif- 
ferent from the local case. We consider non-Gaussianities which effects on the halo clustering is well described by the 
equilateral and enfolded templates. We consider the angular power spectra of these three LSS tracers, as well as their 
cross-correlation power spectra with the CMB temperature fluctuation. Finally, we consider the cubic correction on 
the halo bias, which can be motivated in scenarios like the curvaton model, in which a large cubic correction can be 
produced while simultaneously keeping the /nl correction small p3l . [24| . 

The structure of the paper is as follows. In section [IT] we review the effects of primordial non-Gaussianity on the 
power spectrum Ce- Section |HT1 contains the analysis of power spectrum Ce for NVSS radio sources, SDSS DR6 QSOs 
and MegaZ-LRG. In section ITVl we present the method used to derive constraints on /nl- Section fVl contains our 
main results. We conclude with a discussion in section IVT1 



II. EFFECTS OF PRIMORDIAL NON-GAUSSIANITY 



The effects of non-Gaussianity on the source clustering properties arise because a non-zero /nl affects the halo 
mass function and enhances the halo clustering on large scales. The second effect is the dominant one on large scales. 

In the presence of primordial non-Gaussianity, the mass function uncI-^, z , /nl) can be written in terms of the 
Gaussian one riQ m (M, 2), for which a good fit to the results of simulations is provided by e.g., the Sheth-Tormen 
formula (2f| , multiplied by a non-Gaussian correction factor [24|, [2(| [2?| 2 : 



R NG (M,z, /nl) = 1 
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where the normalized skewness of the density field, S^m, is oc /nl, cm denotes the rms of the dark matter density 
field linearly extrapolated to z = and smoothed on scale R corresponding to the Lagrangian radius of a halo 
of mass M, and <5 CC is the critical density for ellipsoidal collapse, calibrated on N-body simulations (30l - [32l |. For 
high peaks (c^oc/cm 3> 1) and small /nl, S oc is slightly smaller than the critical density for spherical collapse, 
5 c (z) = A c (z)D(0) / D(z) where D(z) is the linear growth factor, and A c (z) ~ 1.68 evolves very weakly with redshift. 
The large-scale halo bias is also modified by the presence of non-Gaussianity [9MTTI l30l . |H, [34[: 



&NG(z) - &g(z) — 2(6g(«) - l)/NL^cc(z)aM(fc) 



(3) 



2 Although attempts have been made to derive directly an expression for the non-Gaussian mass function [28, 29]. 
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where the factor aivi(fc) encloses the scale and halo mass dependence. Here, we consider three types of non-Gaussianity. 

In practice, we find that, on large scales, Q!m(&) c< l/k 2 , is independent of the halo mass for the local type. The 
factor aM(i;) for the enfolded template is proportional to l//c, while the equilateral template gives an almost scale- 
independent ajy(i). It is important to note here that the equilateral and enfolded templates were created to have 
factorizable expressions that gave a good description of the actual bispectra, on average, over all configurations. For 
the the halo bias effect, nearly squeezed configurations dominate. As discussed e.g., in Ref. |IH the templates may not 
necessarily offer a good approximation of the actual bispectra generated by higher-order derivatives or vacuum-state 
modifications. For all inflationary single-field models the bispectrum scales either as l/k or as l/k 3 in the squeezed 
limit (i.e. either like the equilateral template or like the local case) [34[. The orthogonal template instead scales like 
l/k 2 in the squeezed limit. Nevertheless it is interesting to consider such an intermediate case: inflationary models 
that go beyond single field can show a range of scaling in the squeezed limit (e.g., ref. Hf|). 

We assume that the large-scale linear halo bias for the Gaussian case is [251 ]: 
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where z? is the halo formation redshift, and z Q is the halo observation redshift. As we are interested in massive haloes, 
we expect that Zf — z Q . H ere, q = 0.75 and p = 0.3 account for the non-spherical collapse and are a fit to numerical 
simulations (see also refs. |36l438j ). Finally, the weighted effective halo bias is given by 



b^ G (M min7 Z 7 kJ NL ) 



Jm ■ &NG"NGdM 



being the minimum halo mass hosting a source of the kind we are considering. 



(5) 




FIG. 1: Normalized redshift distributions, dN/dz, of the different types of objects considered for our analysis. NVSS radio 
sources (red dashed line), SDSS DR6 QSOs (black solid line), and SDSS DR6 LRGs (blue dash-dotted line). 



III. OBSERVED POWER SPECTRUM 



A. Discrete Sources Maps 



In this subsection we describe three LSS probes, NVSS radio sources, high redshift SDSS DR6 QSOs and McgaZ 
DR7 LRGs sample. In figurc[T]we show their redshift distribution, dN/dz(z). All distributions are normalized to unit 
integral. 
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1. NVSS Radio Sources 

The NVSS [l4| offers the most extensive sky coverage (82% of the sky to a completeness limit of about 3 mJy 
at 1.4 GHz). Wc start by confining our analysis to NVSS sources brighter than 10 mJy, since the surface density 
distribution of fainter sources suffers from declination-dependent fluctuations [39j]. Density gradients in the NVSS 
catalog become increasingly unim por tant as the source brightness threshold is increased. Figure 6 of ref. [39j . figure 
1 of ref. |4(j and figure 3 of ref. [4ll | showed variations of NVSS source surface density as a function of declination 
for sources with different flux thresholds. Dim sources are strongly affected by density gradients, however this does 
not appreciably happen for the brighter sources. Furthermore, in our previous work we used the NVSS sources with 
different flux thresholds S > 10 mJy and S > 20 mJy to compute the ACF of NVSS source and its CCF with CMB 
temperature fluctuations and found that the obtained ACF and CCF results are stable [2l| . Also we mask the stripe 
\b\ < 5°, where the catalog may be substantially affected by Galactic emissions. Here, we have explicitly checked that 
the choices of masked stripes with |6| < 10° and \b\ < 20° have negligible impact on our final results. In order not to 
excessively reduce the sample, in our analysis we consider the NVSS sources brighter than 10 mJy and the masked 
strip with \b\ < 5° to reduce the contamination from Galactic emissions. The NVSS source surface density at this 
threshold is 16.9 deg -2 and the redshift distribution at this flux limit has been recently determined by ref. [421] . Their 
sample, complete to a flux density of 7.2 mJy, comprises 110 sources with Si.4GHz > 10mJy, of which 78 (71%) have 
spectroscopic rcdshifts, 23 have redshift estimates via the K-z relation for radio sources, and 9 were not detected in 
the K band and therefore have a lower limit to z. We adopt here the smooth description of this redshift distribution 
given by ref. [43| , which is shown in the red dashed line of figure [TJ 

dN/dz = 1.29 + 32.37z - 32.89z 2 + 11.13z 3 - 1.25z 4 . (6) 

The mean redshift of this sample is z = 1.23 and the redshift range is < z < 3.5. 



2. SDSS DR6 Quasars 

The SDSS DR6 QSO catalog released by Ref. J15| contains about 10 6 objects with photometric redshifts ranging 
from 0.065 to 6.075 over a total area of 8417 deg (~ 20% of the whole sky). We refer the reader to ref. [l5| for a 
detailed description of the object selection with the non-parametric Bayesian classification kernel density estimator 
(NBC-KDE) algorithm. We use the electronically-published table that contains only objects with the "good" flag 
with values in the range [0,6]. The higher the value, the more probable for the object to be a real QSO (see §4.2 
of ref. [IH for details). Furthermore we restrict ourselves to the "uvxts=l", i.e. to QSOs clearly showing a UV 
excess which should be a signature of a QSO spectrum (in this case we have N qso sa 6 x 10 5 QSOs). As we know, 
several systematics effects, including the galactic extinction by dust, sky brightness, number of point sources and 
poor seeing, could potentially affect both the observed auto-/cross-correlation power spectra. In our previous work we 
carefully checked for their contribution in calculations and in particular we considered extinction and point sources 
contamination. We found that the dominant systematic effect is the extinction 0, |46| . Therefore, in order to 
minimize the effect of Galactic extinction on the observed QSO distribution, we use the extinction mask A g < 0.18 
only in our following analysis f44T - |46j ] . These masks will remove about ~20% of the considered area. 

We fit the redshift distribution dN/dz of the DR6 QSO sample with a function of the form [44[ (black solid line of 
figure [1|): 



dN _ p z m 
-feW ~ ^sj+l) ex P 



(7) 

The best-fit values of the parameters are m = 2.00, j3 = 2.20, zq = 1.62; the mean redshift of the sample is z ~ 1.49. 



3. MegaZ DR7 LRG 

We use the updated MegaZ LRG DR7 sample [22|, which contains ~ 1.5 x 10 6 galaxies from the SDSS DR7 in the 
redshift range 0.4 < z < 0.7 with limiting magnitude i < 20. To reduce stellar contamination, there is a variable of 
the MegaZ neural network estimator <5 sg , defined such that S sg = 1 if the object is a galaxy, and S sg = if it is a 
star. For a conservative analysis, we choose a cut <5 sg > 0.2, which is reported to reduce stellar contamination below 
2% while keeping 99.9% of the galaxies. In addition to the SDSS DR7 geometry mask, we also add two foreground 
masks to account for seeing (removing pixels with median seeing in the red band larger than 1.4 arcsec) and reddening 
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(removing pixels with median extinction in the red band A r > 0.18) [471 ]. The redshift distribution function (blue 
dash-dotted line of figure [1} in this case is found directly from the photometric redshifts that are given in the catalogue 
1481 . which is similar to a Gaussian function: 



dN, , 1 

-T~W = fR — ex P 
dz y/2ir a 



(8) 



2 a 2 

where the mean redshift zq = 0.55 and the deviation a = 0.055. 

B. The Power Spectrum Measurement 

1. The pseudo-Ci estimator 

A two dimensional density field er(n) defined over the full sky can be decomposed in a series of spherical harmonics 
Yim and their corresponding coefficients ai m : 

oo £ 

o-(n) = ^2 E a e™Ytrn(n) , (9) 

t=0 m=-£ 

with 

a lm = J dn n a(n)Y e * m (n) . (10) 

On the full sky spherical harmonics are a complete and orthogonal basis set. However, in practice, a masked (incom- 
plete) sky is observed. Thus the observed di m coefficients on the partially cut sky become (49| : 

at m = f dn n a(n)W(n)Y e * m {n) , (11) 

where VF(n) is the position dependent, weight function imposed by the mask and, optionally, the adopted weighting 
scheme. From these coefficients we can obtain the observed pseudo power spectrum Ci, the pseudo-CV estimator [4Q| : 

1 l 

The pseudo power spectrum Ci, given by the weighted spherical harmonic transform of a masked map, is clearly 
different from the full sky true power spectrum Ce- The expectation value of Ci is related to the true Ci by a 
convolution: 

(C e ) = Y, c t'Gu> , (13) 

v 

where the coupling matrix Gw, describing the mode coupling resulting from the weight function W(n) [5(j, can be 
expressed in terms of 3j symbols as: 

G «, = (24 + D "f "Ji±lm,(y^A 2 . (14, 

e 3 =\ti-tx\ 

We is the angular power spectrum of the weight function: 

1 1 

m = 27TT E l^'"l 2 ' ( 15 ) 

m=l 

where the weight coefficients wi m are: 

wim = J dn n W(n)Y e * m (n) . (16) 
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For the LSS catalogues presented above, in practice, we construct the pixelized maps using the HEALPix software 
package [5l[ with high resolution N s[dc = 512, yielding pixel areas of 6.87' x 6.87' and estimate their pseudo power 
spectrum Ci using the HEALPix function "anaf ast". We also subtract the shot noise contribution ACl/N from the 
observed pseudo power spectrum, where Af2 is the surveyed area and N is the observed number of sources. The 
pseudo power spectrum is measured up to £ max = 400, at these multipoles the magnitude of the estimated error bars 
is dominated by shot-noise. In order to avoid the effects of gauge corrections on the power spectrum on very large 
scales [HI, HH, we set the minimal multipole £ m j n = 10. The effects of non-Gaussianity are large on very large scales. 
By ignoring I < 10 we are neglecting the scales where the signal is maximal. However the theoretical predictions for 



the measured correlations is uncertain on very large scales as discussed in ref. [531 ] . In ref. [21[ we estimated that 
these uncertainties could introduce a systematic error on /nl of 5. Here we take a more conservative approach and 
simply ignore the largest scales, even if this implies a reduction of the signal-to-noisc. 

2. Covariance matrix of the pseudo-C \ estimator 
In the signal dominated limit, the inverse of the power spectrum covariance matrix Fa> is (53 |: 

pmask = (2l+l)G«, 

*»> 2(C t + N e )(C e , + N t >) ' [ ' 

where Ci is the theoretical angular power spectrum, Ng = AO/A denotes the shot noise term. Ci and Ni are 
the necessary error contributions from both cosmic variance and shot noise, respectively. For the full sky survey, 
Gtxi 2 = ^1^2 /sky and the covariance matrix becomes diagonal: 

pfull _ ( 2 ^ + 1 )/sky c /-.ox 

F "' ~ 2(C e + N e )(C,+N e ,) Su ' ■ (18) 

3. Cross- correlation power spectrum 

Besides the angular power spectrum of each LSS tracer, we also consider the cross-correlation power spectrum 
among them. In figure [U we show that the redshift distributions of these three tracers overlap. We can thus expect 
some cross-correlation signal among them. We also use the HEALPix software to estimate the cross-correlation power 
spectrum Cf Y between tracer X and Y: 

^ XY - ^TT E a£.*£ > ( 19 ) 

where a\ m are the observed coefficients of tracer i and they are computed on the common sky area only. Consequently, 
the inverse of the covariance matrix becomes: 

f mask = {2Z+l)G U > ^ 

" ^{Cf Y Y + (Cf X + Af X )(C/ Y + NJ Y ) yj (CfY)2 + (C XX - N XX ){C YY + JVYY) ' 

where the shot noise term N\ = Afi/A; is scale independent, Ail is the common area of two surveys, Ni is the 
observed number of sources in this common area. 

Similarly, we also estimate the cross-correlation power spectrum C* XT between the LSS number density map and 
the CMB map. While the Planck satellite is operating (5f|, here, we use the the 7- years Internal Linear Combination 
(ILC) map of the CMB provided by the Wilkinson Microwave Anisotropy Probe (WMAP) team [5(|. Since an 
estimate of the uncertainty associated to this map is not provided, we have checked that using other CMB maps (e.g. 
the 5 years ILC map by the WMAP team) does not change the results in any significant way. For the cross-check, 
we also use the V and W band all-sky 7-years CMB maps provided by WMAP (531 and smooth these maps to the 
resolution of ILC map (60 arcmins). Using these maps, we perform all the calculations again and find that our results 
are undistinguishable. We adopt the WMAP KQ75 mask, excluding about 30% of the sky at low Galactic latitude, 
to avoid most of the residual Galactic contamination. 

We still use the estimator of cross-correlation power spectrum Cf T : 

1 1 

— 2( ZjZ l ®tm®tm ' (21) 
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where aj m are the observed coefficients of the CMB map AT(n), also computed on the common sky area only. The 
inverse of the covariance matrix is: 

pnaek _ (21 + \)Gw ^ 

" sJiCT) 2 + (Cf x + N™)(CJ T + iV7 T )^(C* T )2 + (C$ x + N™)(Cjr + N™) ' 

where Cj T is the theoretical angular power spectrum of CMB map, Nj" 1 is measurement error of the CMB temper- 
ature power spectrum, in which the main source is the pixel noise which accounts for the imperfections in measured 
temperature on the sky. Pixel noise for WMAP is reported as er p i X = <tq/ \/Nob B , where co is the noise per observation 
and iVobs is the number of observations per given pixel [58$ . After getting the theoretical Cj and C XT , we include the 
correction of Gaussian beam B(£) (using the HEALpix gaussbeam routine): Cj T = Cj T B(£) 2 , Cf T = Cf T B(£). We 
neglect the pixel window function since its contribution at A S id c = 512 is negligible at the resolution of the WMAP 
ILC map which is 60 arcmins FWHM. 

IV. METHOD & DATA ANALYSIS 

We use the publicly available package CAMB_sources 3 [59[ to calculate the theoretical power spectrum. In our 
analysis, we use the Halofit [60] built-in routine for non-linear correction to obtain the fully-evolved, nonlinear matter 
power spectrum P(k) at any epoch. We perform a global fitting of cosmological parameters, including /nl, for the 
data of section Hill including also the datasets described below, using the CosmoMC package [6l|, a Markov Chain 
Monte Carlo (MCMC) code, modified to calculate the theoretical ACF and CCF. We assume purely adiabatic initial 
conditions and a flat Universe, with no tensor contribution to primordial fluctuations. The following six cosmological 
parameters are allowed to vary with top-hat priors: the dark matter energy density parameter Sl c h 2 £ [0.01,0.99], 
the baryon energy density parameter Fl^h 2 £ [0.005,0.1], the primordial spectral index n s £ [0.5, 1.5], the primordial 
amplitude log[10 10 A s ] £ [2.7,4.0], the ratio (multiplied by 100) of the sound horizon at decoupling to the angular 
diameter distance to the last scattering surface 9 S £ [0.5, 10], and the optical depth to reionization r £ [0.01,0.8]. 
The pivot scale is set at fc s o = 0.05 Mpc -1 and do not consider massive neutrinos and dynamical dark energy. Besides 
these six basic cosmological parameters, we have two more parameters related to the power spectrum Ct data: the 
non-Gaussianity parameter, the minimal halo mass Af m j n . When we combine all the different probes together we have 
four parameters: three minimal halo masses for three LSS surveys and the non-Gaussianity parameter. 

The model powe r sp ectrum Cf 1 is compared with the observed values CV, respectively, through the Gaussian 
likelihood function [62 ]: 

ln£ cx -\Y^ G ^ ~ C t )F w {C$ ~ C v ) , (23) 
w 

where the Fisher matrix Fui is the inverse of the power spectrum covariance matrix. While it is true that at very low 
I a Gaussian likelihood is not a good approximation for CMB data, at high I the Gaussian approximation works well 
even for CMB, for LSS it is customary to adopt a Gaussian likelihood. 

The following cosmological data are also included in the fit: i) power spectra of CMB temperature and polarization 
anisotropics; ii) baryonic acoustic oscillations (BAOs) in the galaxy power spectra; iii) SNIa distance moduli. 

To deal with the 7-years WMAP (WMAP7) CMB temperature and polarization power spectra we use the routines 
for computing the likelihood supplied by the WMAP team [63j . The WMAP7 data are used only to improve the 
constraints on the six basic cosmological parameters, not to constrain /nl- 

The BAOs |64| can, in principle, measure not only the angular diameter distance, Da{z), but also the expansion 
rate of the Universe, H(z). However, the limited accuracy of current data only allows us to determine the ratio 
between the distance scale defined by ref. [65| : 



D v (z) 



2 n 2 ' * CZ 



"I 1/3 



(1 + 



H(z) 



(24) 



and the comoving sound horizon at the baryon-drag epoch, r s (zd) (see ref. |66j|). Accurate determinations of the 
distance ratio r s (zd)/ D v (z), r s {zd) have been obtained by Ref. [641 ] : 

r s (z d )/D v (z = 0.20) = 0.1905 ±0.0061, 



Available at http://camb.info/sourccs/ 
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r s (z d )/D v (z = 0.35) 



0.1097 ±0.0036. 



(25) 



We adopt these values as Gaussian priors. 

The SNIa distance moduli provide the luminosity distance as a function of redshift D^(z) which provides strong 
constraints on the dark energy evolution. In this paper we will use the latest SNIa data sets from the Supernova 
Cosmology Project, "Union2 Compilation" which consists of 557 samples and spans the redshift range < z < 1.55 
I67I . In this data set, they improved the data analysis method by using and refining the approach of their previous 
work [68j . When comparing with the previous "Union Compilation" , they extended the sample with the supernovae 
from refs. H [6§Ei. The authors also provide the covariance matrix of data with and without systematic errors 
and, in order to be conservative, we include systematic errors in our calculations. In the calculation of the likelihood 
from SNIa we marginalize over the nuisance parameter as done in refs. [721.1731]. 

Furthermore, we add a prior on the Hubble constant, Hq = 74.2 ± 3.6 km/s/Mpc given by ref. [741 . Finally, in the 
analyses we set the minimal halo mass at M min > 10 12 /i -1 M Q consistent with the results of ref. [75j . 
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FIG. 2: Observed angular power spectra of three LSS tracers. The black dashed lines are the best fit models using angular 
power spectra only. The red solid lines are the best fit models /nl = 48 when using all data together. For illustration purposes, 
we show the binned power spectra with the bin size A£ — 10. 



V. NUMERICAL RESULTS 



In this section, we present constraints on the primordial non-Gaussianity parameter /nl from the observed LSS 
power spectra, in combination with the external data sets introduced above: the WMAP7, BAO and SNIa. These 
extra datascts are only used to constrain the underlying cosmology. In table U and figure 2] we show the constraints 
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FIG. 3: Left Panels: Observed cross-correlation power spectra among three LSS tracers. Right Panels: Observed cross- 
correlation power spectra between CMB ILC map and three LSS tracers. The red solid lines are the best fit models /nl = 48 
when using all data together. For illustration purposes, we show the binned power spectra with the bin size A£ = 10. 



on the primordial non-Gaussianity parameter from different data combinations, after marginalizing over all the other 
parameters. We also plot the observed power spectra data £Ci in figurc[5]and figurc[31 together with the corresponding 
theoretical best fit models for comparison. For illustration purposes, here we bin the observed power spectra with the 
bin size A£ =10. 
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A. Local Type 

1. Angular Power Spectrum 

We begin by considering the constraint on the local type of non-Gaussianity /nl from the angular power spectrum 
only. In the first panel of figure [5J we plot the observed NVSS angular power spectrum data, which is consistent 
with the previous work [76|. When using the NVSS Ct data only, we obtain the marginalized constraint on the 
non-Gaussianity parameter: /nl = 78 ± 52 (ler) and —34 < /nl < 187 (2 a), which is consistent with zero. The 
minimal halo mass, M m j n = 10 12 ' 44±0 26 /i _1 M Q (1 a), turns out to be remarkably close to our previous work [2(j- We 
show the best fit model in this panel (black dashed line) . The amplitude of the power spectrum on large scales (small 
£) is enhanced to fit the data points, by the positive value of /nl- In our previous works 0, [2l|, we used the auto- 
correlation function (ACF) of NVSS radio sources to constrain /nl and found that the NVSS ACF yields a positive 
/nl at more than 2 a. The reason for this difference is that we only use the power spectrum data at £ > 10, and the 
power spectrum approach does not include the information on the total number of NVSS radio sources A to t (which 
was folded in the ACF analysis). If we add a prior on the total number of sources N tot : 0.5 < iVt pt /N ?^ ss < 1.5, we 
obtain a limit of /nl = 62 ± 30 at 1 a. This result is perfectly consistent with previous works |20l [2l|. 

Next, we use the estimated angular power spectrum of the SDSS DR6 QSO sample, which is shown in the second 
panel of figure [2J There is an obvious excess power at large scales, which favors a positive /nl: /nl = 62 ± 26 (1 a) 
and 5 < /nl < 115 (2c). Similarly to our previous work (2lj . we also find evidence for positive /nl at more than 
2 a. However, ref. [18| . using SDSS QSO data, obtained non-Gaussianity constraints consistent with zero at 95% 
confidence level. As we discussed before, ref. [18( used an extension of the DR3 QSO sample that include sources 
that subsequently were released with the DR6 sample: we use the more complete and better calibrated final official 
SDSS DR6 QSO catalog release [U]. Ref. [3] found that their QSO sample at z < 1.45 (QSO0) seem to suffer from 
contamination which they attribute to systematic calibration errors, discarded this sample from their analysis and 
only used the high redshift sample (QSOl). We have checked for this effect and find no evidence in the sample we 
use that the z < 1.45 QSOs have different contamination than the z > 1.45 [2l|. In addition, we remove three narrow 
stripes in the southern hemisphere of the galaxy and set the minimal multipole £ m m = 10, which could reduce the 
effect of star contamination further. 

Finally, we consider the MegaZ DR7 LRG sample. Ref. [22j split this sample in four redshift bins according to 
each source photometric redshift estimate and found a large excess of power over the lowest multipoles {I < 20) in 
the angular power spectrum, which is also found in their previous analysis using MegaZ DR4 LRG sample [771 ]. In 
our analysis, although we do not split the LRG sample in redshift bins, we still find this excess power on large scales, 
which is shown in the third panel of figure [5] There are several explanations for this excess power that range from 
systematic errors as due to contaminants to new physics [78j]. Ref. [78| checked some possible systematic errors of 
this MegaZ LRG sample in details and found that systematics errors do not seem to be responsible for the excess 
power at large scales. We also perform a similar check and find no significant changes in the power spectrum at large 
scales. One of the possibilities is that this excess power could be induced by the primordial non-Gaussianity. Here, we 
use the observed Cg data of MegaZ DR7 LRG sample to constrain the primordial non-Gaussianity for the first time. 
We calculate the angular power spectrum on largest scales and account for the redshift distortion power as described 
in Ref. (79j ). although we set the minimal multipole l m i-a = 10. After marginalizing over other free parameters, we 
obtain the constraint: 

/ NL = 165±105 (la) , -81< / NL < 351 (2a) , (26) 

which is consistent with zero at 2 a. We also obtain a limit of M min = W 12A5±a - 07 h~ 1 M Q at 68% confidence level. 
Using cq. wc could obtain the limit of effective halo bias at the mean redshift of the survey: b(z = 0.55) = 
1.93 ± 0.06 (1 a), which is consistent with other works [8(|, although in our analysis we directly use the full sample. 

When we combine these three observed angular power spectra data together to constrain non-Gaussianity, we obtain 
a limit of /nl = 68 ± 22 at 68% confidence level and 22 < /nl < 108 at 95% confidence level, which implies the 
current angular power spectra of LSS tracers favor /nl > at the ~ 3 a level, since non-Gaussianity adds power on 
large angular scales yielding a good fit to the observed data points. 

As can be seen from figure [2] and table HI there is not a specific measurement that is driving the signal: all the 
different probes yield very consistent results and when combined the error bar on /nl is reduced to yield a > 2a 
result. 
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2. Adding Cross-Correlation Power Spectrum 

Comparing with the angular power spectrum, the constraining power of cross-correlation power spectrum on /nl 
is much weaker [12]. In the panels of figure [31 we plot the observed cross-correlation power spectra Cf Y among three 
LSS tracers and Cf~ T between the CMB map and three LSS tracers. These cross-correlation data points are consistent 
with zero within 1 a error bar but the best fit value /nl = 48 is an excellent fit to the data points, as shown in figure 
[3J However, cross-correlation data give a larger error bar for /nl than that from the angular power spectra data and 
are thus much less constraining. 

First, we consider the single tracer case in combination with its cross-correlation with the CMB map. When we 
combine the angular power spectrum and its cross-correlation power spectrum Cf~ T with CMB map to constrain the 
primordial non-Gaussianity, we obtain /nl = 74 ± 40, /nl = 59 ± 21 and /nl = 153 ± 95 at 68% confidence level 
for NVSS radio sources, SDSS DR6 QSOs and McgaZ DR7 LRGs, respectively. The results are consistent with those 
obtained from the angular power spectrum only, but those cross correlations does not seem to add much information. 

If we consider all the angular and cross-correlation power spectra data sets of three LSS tracers together 4 , we 
obtain the constraint on the local type primordial non-Gaussianity: 

/nl = 48± 20 (la) , 5 < /nl < 84 (2a) . (27) 

This result is compatible with previous estimates fLsl IM I8H486] . and with the WMAP7 limits [63|. The LSS data still 
implies a positive /nl at more than 2 a. In figurc[2]and figure|3j we plot the theoretical best fit model /nl = 48 of our 
non-Gaussianity calculations (red solid lines). These curves match the observed power spectra very well, especially 
on large scales. 

B. Other Shapes 

Generally, the correction to the standard halo bias due to the presence of primordial non-Gaussianity is Abh/bh = 
5 c f3(k), where the function (3(k) is given by [87| : 

W = 8 2 m I dkihMnih) f d»M R (Va) B * {k *'tf*' k) , (28) 

87T ! CT£ I .MR.(fc) J J -I P${k) 

where .Mr = 47rp m i? 3 /3 is the halo mass which is related to the scale R, P$(fc) denotes the expression for the 
primordial bispectrum of the Bardeen potential Pj> its power spectrum and a = k\ + k 2 + 2k\kjJL. 

In the local non-Gaussian case, the expression for the correction to the standard gaussian halo bias simplify [eq. 
(J3J)]. In this subsection, we explore other two different types of primordial non-Gaussianity given by the Equilateral 
and Enfolded templates. 

The equilateral type of non-Gaussianity [H, H2], which can be generated from the inflationary models with 
higher-derivative operators of the inflaton, can be well described by the following template [90]: £>$(fci, fc 2 , fc 3 ) = 
6/ N q L ^ oq (fci,fc 2 ,fc 3 ), where 

F^(k u k 2 ,k 3 ) = (-P#(fci)P*(fca) + 2cyc) - 2[P < ,(k 1 )P^(k 2 )P^(k 3 )} 2 ^ + (P^h) 1 ^ 3 P^(k 2 ) 2 ^ 3 P^(k 3 ) + 5cyc) . (29) 

Then we find that the correction /3(k) in the Equilateral non-Gaussian case is almost scale-independent (see figure 1 
of ref. [13 )• In principle there is a mass dependence (see discussion in ref. [Hj]), but the template and the actual 
physical bispectrum yield different effective /nl (denoted by /nl) for the halo bias. Here, therefore, we consider a 
phcnomcnological scale-independent scaling and neglect the mass dependence. The interpretation of the resulting 
constraint should be interpreted as a constraint on an effective /nl, its connection to the physical bispectrum depends 
on the detailed model of inflation under consideration. 

Using all datasets, we obtain the constraint on the Equilateral type of primordial non-Gaussianity: 

/^ q L = 50 ± 265 (la), - 419 < /£ q L < 625 (2 a) , (30) 



4 In our analysis, we do not take into account the possible correlation between the auto-correlation and cross-correlation power spectra. 
Neglecting this correlation should have a negligible effect. In fact when the cross-correlation signal is much weaker than the auto- 
correlation signal, as it happens here, their covariance is diagonal-dominated with negligible off-diagonal terms. This is well known to 
happen for example when considering the ISW effect from the galaxy-CMB cross correlation and the CMB temperature 
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which is consistent with the WMAP7 results 63 1 



The template of another type of non-Gaussianity "Enfolded" , which arises in models with non-Bunch-Davies initial 
state [H, |9l| or in effective field theories of inflation with higher-derivative interactions [§2|,[93| (see also ref. [§I|), is 
given by ref. |H: B*(ki,k 2 ,k 3 ) = 6f^F cnl (k u k 2 , fc 3 ), where 

F cnf (fc l5 fc 2 ,fc 3 ) = (P*(fci)P*(fc 2 ) + 2cyc) + 3[P$(fc 1 )P <I ,(fc 2 )P 4 ,(fc3)] 2/3 - (P$(fci) 1/3 P*(fc 2 ) 2 / 3 P 4 ,(fc 3 ) + 5 C 2/ C ) . (31) 

In figure 1 of ref. |87[, it can be seen that the correction in this type is proportion to on large scales. 

No single field inflationary model has this scaling in the squeezed limit: all inflationary single field models considered 



so far have only the local or the template equilateral scaling in this limit [341 ]. This is nevertheless an interesting 
intermediate case to consider. In fact non-single field models can yield all the intermediate scalings in the squeezed 
limit, see e.g., ref. J35|. We therefore consider a phenomenological scaling like which could cover this Enfolded 
template case and neglect the mass dependence. In this case, the whole data sets give the constraint on the effective 
non-Gaussianity parameter for the enfolded template: 

/S n L f = 183 ± 95 (1 a) , - 12 < < 358 (2 a) . (32) 

In the WMAP7 paper (63|, a constraint on the Orthogonal type of non-Gaussianity (8j| is given. The correction 
to the halo bias of this template is also proportional to k^ 1 on large scales, but have the opposite sign and a 
normalization factor of f^ h = — 2/nl with a weak mass dependence which we neglect here. Therefore, based on 
the result of Enfolded template, in our simple analysis, we obtain the constraint on the non-Gaussianity parameter 
for the orthogonal template: /£ r L th = -92 ± 47 (la) and -179 < / N £ h < 6 (2a), The orthogonal type of non- 
Gaussianity in the squeezed limit however scales like the equilateral one and, as shown in ref. [34| its effect on the 
halo bias is almost degenerate with the equilateral non-Gaussianity. The enfolded type of non-Gaussianity ari sing 
for example from modified initial state, in the squeezed limit scales like the local template and as shown in ref. [341 ] 
its effect on the halo bias is almost degenerate with the local non-Gaussianity but with a different normalization 
jmod.m. state ^ jyg^iocai This re-scaling can be used to re-interpret the local constraints into constraints for the 
modified initial state non-Gaussianity. 



C. Cubic Correction <?nl 



Finally, we consider the effect of the cubic correction on the halo bias model. In this case, the Bardeen potential $ 
can be conveniently parameterized up to third order by: 

* = 4> + /nl(0 2 - (<t> 2 )) + W 3 , (33) 

where <?nl is dimensionless, phenomenological parameter. Ref. [96| explored the effect of a cubic correction gNh<f> 3 on 
the mass function and halo bias model in detail and gave the expression of the scale-dependent bias correction when 
/nl = and .g N L ^ 0: 

A6(fc, 5NL ) = iekflNiACOSgjA&Cfc, /nl = 1) , (34) 

where the normalized skewness of the density field = S3,m//nl and the parameter ~ 0.6 obtained by the 

fitting in ref. (96J. Then we obtain the constraint on ^nl from all the data sets: 

,9nl = (5.7 ±3.0) x 10 5 (la), - 1.2 x 10 5 < g NL < 11.3 x 10 5 (2 a). (35) 

The result is compatible with that obtained in refs. (96l. [97|. 



VI. SUMMARY & CONCLUSIONS 



Searches for primordial non-Gaussianity can provide key information on the physical mechanism that generate the 
primordial perturbations and can constrain specific inflationary models. A relatively new and promising method 
to detect primordial non-Gaussianity exploits the fact that the clustering of dark matter halos (which host massive 
galaxies and QSOs) is modified on large-scales in a scale-dependent way. This effect is called non-Gaussian halo bias. 
To access these large scales, large-scale structure surveys covering large volumes, high redshift and tracers of massive 
halos arc most suited. Among the currently available surveys, the NVSS and the SDSS QSO (in particular the DR6 
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TABLE I: 1, 2a constraints on the primordial non-Gaussianity from different data combinations. We report the mean values 
and the Bayesian central credible interval, marginalized over all other parameters. 



Datasets 


Non-Gaussianity 


WMAP7+BAO+SN 


la 


2(7 


Local Type /nl 


+Cf x (NVSS) 
+C| cx + Cf r (NVSS) 


78 ±52 
74 ±40 


[-34, 187] 
[-16,166] 


±Cf x (QSO) 
+Q xx + C] CT (QSO) 


62 ±26 
59 ±21 


[5, 115] 
[17, 103] 


±Cf x (LRG) 
+C XX + C? T (LRG) 


165 ± 105 
153 ± 95 


[-81,351] 
[-51, 347] 


±£7]™ (ALL) 
±Cf x + Cf Y + C XT (ALL) 


68 ±22 
48 ±20 


[22, 108] 
[5, 84] 


Equilateral Template /nl 


+Cf x + Cf Y ± C XT (ALL) 


50 ± 265 


[-419,625] 


Enfolded Template /nl 


+Cf x + Cf Y + C XT (ALL) 


183 ± 95 


[-12,358] 


Cubic Correction jnl X 10 


+Cf x + Cf Y + Cf T (ALL) 


5.7 ±3.0 


[-1.2,11.3] 



NVSSC, - 

NVSS 

c xx + c j 



QSOC 



XX 



QSO 

c xx + c XT 

LRG C, 3 ™ 
LRG 

C™ + c, XT 
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ALL C ( xx + 
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FIG. 4: Median values (red points) and la and 2a limits on the local type of primordial non-Gaussianity /nl obtained from 
different data combinations. The vertical line denotes the Gaussian case /nl = 0. 



catalog) arc the largest. Previous work had detected excess large-scale power in their auto-correlation function yielding 
a detection of a positive value of the non-Gaussianity parameter /nl at > 2a level. Here we have revisited the issue 
using the angular power spectrum rather than the correlation function. While the two statistics should enclose the same 
physical information they are affected by systematics effects in different ways: this complementary analysis offers thus 
a consistency check and a test for possible systematics. We have also used the recently updated MegaZ-LRG survey: 
a previous release of the same survey had also been found to have excess large-scale power, but mostly at £ < 10. In 
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an effort to be as conservative as possible in our analysis, we excuded all multipoles at I < 10: these multipoles are 
affected by theoretical uncertainties and systematics such as calibration, sample contamination etc. These also probe 
the largest scales, thus these are the multipoles where most of the non-Gaussian signal is likely to be concentrated. 
We combine the survey auto power spectra with the cross-power spectra among the surveys and of each of the survey 
with the CMB temperature map. In our analysis we marginalize over the LCDM cosmological parameters with a 
prior imposed by WMAP7 angular power spectra, BAO, supernovae and Hubble constant measurement. We have 
investigated the constraints on the parameter /nl, characterizing primordial non-Gaussianity not only for the popular 
local shape but also for the equilateral, enfolded and orthogonal templates. We have also considered constraints on 
a cubic non-Gaussianity correction (parameterized by the <7nl parameter). We find no evidence for non-Gaussianity 
for any "shape" which halo bias effect on large scale scales less steeply than 1/fc 2 and no evidence for non-zero (?nl- 
However for the local type of non-Gaussianity we find a ~ 2a signal for positive /nl when all data-sets are combined. 
The results are summarized in figure 2J no single data set drives the signal and the central /nl values recovered 
are fully in agreement with previous analyses. Due to the more conservative approach taken here (and thus larger 
error- bars), the only single survey that yield more than 2a detection is the SDSS DR7 QSO survey. The NVSS survey 
central recovered /nl value is also in agreement with both the QSO measurement and previous analyses but the 
error-bar is larger. The combination of the data sets yields: /nl = 48 ± 20 and 5 < /nl < 84 at 1 and 2 a. Thus 
the tantalizing hint for a positive local /nl found from the auto-correlation function remains despite the different 
analysis. We have investigated several possible sources of systematic errors such as the sources number density and 
the contamination by Galactic emissions for NVSS sources; different choices for the CMB temperature fluctuations 
templates; possible contamination of stars for the SDSS and LRG samples. Our main results are stable when these 
systematic errors are considered. Such a result would have profound implications for inflationary mechanisms. The 
improved statistical power of forthcoming Largc-Scalc Structure surveys is needed to fully resolve the issue. Moreover, 
it will be also important to use the CMB data which will be provided by Planck. 
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